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Abstract 

The composition and equation of state of the outer crust of non-accreting neutron stars is 
computed using accurate nuclear mass tables. The main goal of the present study is to understand 
the impact of the symmetry energy on the structure of the outer crust. First, a simple "toy model" 
is developed to illustrate the competition between the electronic density and the symmetry energy. 
Then, realistic mass tables are used to show that models with a stiff symmetry energy — those 
that generate large neutron skins for heavy nuclei — predict a sequence of nuclei that are more 
neutron-rich than their softer counterparts. This result may be phrased in the form of a correlation: 
the larger the neutron skin of ^'^'^Pb, the more exotic the composition of the outer crust. 

PACS numbers: 26.60.Kp, 26.60.Gj, 21.65.Ef, 21.10.Dr 
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I. INTRODUCTION 



Neutron stars are gold mines for the study of nuclear systems under extreme conditions of 
density and isospin asymmetry Spanning many orders of magnitude in density, neutron 

stars display exotic phases that cannot be realized under normal laboratory conditions. 
While the most common perception of a neutron star is that of a uniform mantle of neutrons 
packed to densities that may exceed that of normal nuclei by up to an order of magnitude, 
the reality is far different and much more interesting. First, although the uniform liquid 
mantle (also known as the outer core) is indeed composed mostly of neutrons, a small 
fraction of protons and an equal number of charged leptons {i.e., electrons and perhaps 
even muons) must be present to maintain beta equilibrium. The precise proton fraction in 
the neutron star is controlled by the symmetry energy, a quantity that imposes a penalty 
on the system as it departs from the isospin symmetric limit of equal number of neutrons 
and protons. Second, at densities that are below nuclear matter saturation density the 
uniform phase becomes unstable against density fluctuations. This non-uniform region of 
the neutron star constitutes the crust, which itself is divided into an inner and an outer 
region (see Fig. [1]). In the outer crust — the main focus of the present study — the system 
is organized into a Coulomb lattice of neutron-rich nuclei embedded in a uniform electron 
gas [31]. As the density increases, nuclei become progressively more neutron rich until the 
neutron drip region is reached; this region defines the boundary between the outer and the 
inner crust. As in the case of the outer crust, the inner crust also consists of a Coulomb 
lattice of neutron-rich nuclei embedded in a uniform electron gas. Now, however, a uniform 
neutron vapor permeates the system. As the density continues to increase in the inner 
crust, the system is speculated to morph into a variety of complex and exotic structures, 
such as spheres, cylinders, rods, plates, etc. — collectively known as nuclear pasta 0, 0]. 
As the density increases even further, uniformity is eventually restored at about one third 
of normal nuclear matter saturation density. Finally, at ultra high densities it has been 
established that the ground state of hadronic matter becomes a color superconductor in a 
color-flavor-locked (CFL) phase jl, 0]- It is unknown, however, if the density at the core of 
a neutron star may reach the extreme values required for the CFL phase to develop. Thus, 
other exotic phases — such as meson condensates, hyperonic matter, and/or quark matter 
— may be more likely to harbor the core of neutron stars. Figure [T] is believed to represent 
a plausible rendition of the structure of a neutron star. 

As stated earlier, the main focus of our present study is the outer crust of the neutron 
star. In particular, we are interested in studying the sensitivity of the composition of the 
outer crust to the model dependence of the symmetry energy. The outer crust comprises a 
region spanning about seven orders of magnitude in density; from about lO^g/cm^ up to a 
neutron-drip density of about 4 x lO^^g/cm^ [3]. Although small relative to nuclear matter 
saturation density (2.5 x lO^^g/cm^), at these densities the electrons (present to maintain 
charge neutrality) are no longer bound to nuclei and move freely throughout the crust. 
Moreover, at these low nuclear densities it is energetically favorable for the nuclei to arrange 
themselves in a crystalline lattice. At the lowest densities, the electronic contribution is neg- 
ligible so the Coulomb lattice is populated by ^^Fe nuclei. However, as the density increases 
and the electronic contribution becomes important, ^^Fe ceases to be the most energetically 
favorable nucleus. Instead, it becomes energetically advantageous for the system to lower 
its electron fraction by having the energetic electrons capture onto protons, with the excess 
energy carried away by neutrinos. The resulting nuclear lattice is now formed by nuclei 
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having a slightly lower proton fraction than ^^Fe {e.g., ^^Ni). As the density continues to 
increase, the nuclear system evolves into a Coulomb lattice of progressively more neutron- 
rich nuclei until the critical neutron-drip density is reached. The essential physics of the 
outer crust is then nicely captured by a "tug-of-war" between an electronic contribution and 
the nuclear symmetry energy, with the former favoring neutron-rich nuclei while the latter 
favoring fairly symmetric ones. The neutron-rich nuclei that populate the Coulomb lattice in 
the outer crust are on average more dilute than their more symmetric counterparts because 
of the development of a neutron skin. As a result, these nuclei may become sensitive to the 
symmetry energy below nuclear matter saturation density. However, whereas the symmetry 
energy is relatively well known around saturation density, its density dependence [e.g., its 
slope) is poorly constrained. This may affect the composition of the crust. 

While some theoretical constraints are starting to emerge [sl, [§| , the density dependence 
of the symmetry energy remains largely unknown. Indicative of this fact is that accurately- 
calibrated models of nuclear structure (both relativistic and non-relativistic) that reproduce 
a variety of ground-state properties across the periodic table differ significantly in their 
predictions for the density dependence of the symmetry energy. Yet these same models have 
been used to uncover a strong correlation between the pressure of pure neutron matter at 
saturation density and the neutron skin of heavy nuclei: the larger the pressure the larger 
the neutron skin 0, [HI. (Note that the pressure of pure neutron matter equals that of 



the symmetry energy at saturation density.) This fact may be illustrated using the two 
accurately-calibrated models that will be employed throughout this manuscript, namely, 
NL3 [H, and FSUGold Whereas NL3 predicts a pressure at saturation density of 
Po~6 MeV/fm^ and a corresponding neutron skin in ^osp]-, Rn — Rp^O.28 fm, FSUGold 
predicts the significantly lower values of Po~3 MeV/fm and Rn—Rp f^0.21 fm, respectively. 
The upcoming Parity Radius Experiment (PREx) at the Jefferson Laboratory will provide 
a unique experimental constraint on the density dependence of the symmetry by measuring 
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the skin thickness of ^'^^Pb accurately and model independently via parity- violating electron 
scattering [ 
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The correlation between the density dependence of the symmetry energy 
and the neutron skin of heavy nuclei opened new horizons in nuclear astrophysics. Novel 
correlations between the neutron skin of ^osp]^ g^^^ myriad of neutron-star observables were 
developed as a result of the similar composition of the neutron skin of a heavy nucleus and 
the inner-crust /outer-core of a neutron star 17, 18, 2^, 21 1. One particularly interesting 
correlation of direct relevance to the crustal region is a "data-to-data" relation between the 
neutron skin of ^°^Pb and the crust-to-core transition density [l7j . 

The recent observation of intense pulses of energetic gamma rays followed by fainter 
periodic signals emitted from highly magnetized neutron stars (or "ma^netars") are sure 



to provide an additional new tool in the study of neutron-star structure [22|, |23|, |2J, [25 
The discovery of high-frequency oscillations in the tails of giant flares from Soft Gamma 
magnetars with magnetic fields in excess of 10^"^ 



27, 28 
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Repeaters, i.e., magnetars witn magnetic neias m excess ot iU"' gauss 
believed to be associated with seismic vibrations of the neutron star crust. Early theoretical 
models that assume a liquid-core/solid-crust interface suggest torsional shear oscillations 
of the crust as the most likely modes of excitation in a magnetar. The shear modulus of 
the crust acts as a restoring force for these modes and such a structural property is highly 
sensitive to the composition of the crust and, thereby, to the nuclear matter equation of 
state. Indeed, the shear-mode oscillations depend strongly on the neutron star mass, radius. 



and crustal composition — all properties sensitive to equation of state [22(]. Moreover, ratios 



of frequencies with different nodal structures may be used to determine the thickness of the 
crust, an observable highly sensitive to the equation of state and particularly to the density 
dependence of the symmetry energy Hill. Hence, as techniques continue to improve, we 
expect that neutron-star seismology will provide stringent limits on the equation of state of 
neutron-rich matter. 

The manuscript has been organized as follows. The formalism required to compute the 
composition and equation of state of the outer crust is developed in Sec. [Tll In Sec. IIIII we 
employ several realistic nuclear-mass models to compute the structure of the outer crust. 
While not nearly as comprehensive as the recent study performed by Ruester, Hempel, and 
Schaffner-Bielich jlof, ours include a simple "toy model" that provides critical insights into 
the role played by the symmetry energy. Moreover, in the same section we illustrate the 
impact of the density dependence of the symmetry energy on the sequence of neutron-rich 
nuclei present in the outer crust. Our results and conclusions are summarized in Sec. IIV[ 



II. FORMALISM 



In this section we develop the formalism necessary to understand the composition and 
equation of state of the outer crust of a neutron star. The formalism follows closely the 
seminal ideas introduced by Baym, Pethick, and Sutherland back in 1971 3|. For more 
recent references that employ modern nuclear mass tables see Refs. 0, l30l . l3l[ | . The central 
question that one aims to answer is the following: what is the ground state of cold, fully- 
catalyzed matter for densities between complete ionization (p ~ lO^g/cm^) and "neutron 
drip" (p lO^^g/cm^)? Since at these densities uniform matter is unstable against cluster 
formation, a Coulomb lattice of nuclei embedded in a uniform free Fermi gas of electrons 
is formed. Thus, the composition of the outer crust is determined by that nucleus (with 
neutron number A^, proton number Z, and baryon number A = N + Z) that minimizes 
— for each density — the total energy per nucleon of the system. In the outer crust {i.e., 
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before neutron drip) the energy per nucleon consists of three different contributions: nuclear, 
electronic, and lattice. That is. 



e{A,Z;n) = en + ee + Ei , (1) 

where the baryon density is denoted by n = Atotai/V . The nuclear contribution to the total 
energy per nucleon is simple and independent of the density. It is given by 

en{N, Z) = , with M{N, Z) = Nrrin + Zrrip - B{N, Z) . (2) 

Here M{N,Z) is the nuclear mass, B{N,Z) is the corresponding binding energy, and m„ 
and rrip are neutron and proton masses, respectively. 

The electronic contribution — at least at the densities of interest (p > 10''g/cm'^) — is 
modeled as a degenerate free Fermi gas [3]. That is. 
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ee{A,Z;n) = — = — -/ p^/p^ + mldp, (3) 

where S'e, rUe, and ppe are the electronic energy density, mass, and Fermi momentum, re- 
spectively. Note that the electronic Fermi momentum and baryon density are related as 
follows: 

ppe = (Svr^ne) = {37r'yn) ^ y^l'p^ , (4) 

where the electron fraction y = Z/A has been defined. Moreover, for future convenience the 
following definition of the overall Fermi momentum has been introduced: 

PF=(3vr2n)^/^ (5) 

As the integral in Eq. ([3]) may be evaluated analytically, the electronic contribution may be 
computed in closed form. That is, 

4 

£e(A, Z] n) = xfVf [xl + ylj - ln{xF + Vf) , (6) 

where dimensionless Fermi momentum and energy have been defined as follows: 



xf = ^— and yF = — = \/l + x'i. (7) 

nie me V 

We now discuss the last term in Eq. ([1]). Whereas the Coulomb repulsion within the in- 
dividual nuclei has been properly included in Eq. ([2]), the Coulomb repulsion among nuclei 
as well as their interactions with the uniform electron background has not. At the densi- 
ties/temperatures of relevance to the outer crust, namely, large enough for full ionization 
but small enough for the Coulomb repulsion among nuclei to dominate over their kinetic 
energy — Wigner has shown (in the context of the electron gas) that the system will crys- 
tallize into a body-centered-cubic lattice [H, 33, The last term in Eq. ([T]) represents 



the lattice contribution to the energy per particle. The calculation of the potential energy 
of the Coulomb lattice is complicated. It consists of divergent contributions that must be 
canceled as required by the overall charge neutrahty of the system. Fortunately, accurate 
numerical calculations for the electron gas have been available for a long time 35, 3^ and 
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these results can be readily generalized to the present case jl]. Indeed, the lattice energy 
per nucleus may be written as follows: 

^ = -(1.81962)^^ = -(1.79186)^ . (8) 

Nc a 2ro 

where is the number of nuclei {i.e., A-body clusters), a is the lattice constant, and tq is a 
length scale related to the volume per nuclei. For the particular case of a body-centered-cubic 
lattice, these quantities are related in the following way: 

n . f a\ /7r\i/3 

n^a^ = -^a^ = 2 or — = - . (9) 



V \2roJ \3. 

Using the fact that the number of nuclei is related to the total baryon number of the system 
as 

_ ^total _ -^total _ Ztotal H 

the lattice contribution to the energy per baryon may be written in closed form as follows: 

^ , (1.79186) (Ze)2 . 

^M^;n) = -4^^. (11) 

Note that here -Rq refers to the length scale associated to the volume per baryon (ro = 
A^/^Ro). That is. 

Using these definitions, the lattice contribution becomes equal to 

EfXA, Z; n) = -Ce—:rpp (with d = 3.40665 x 10"^) . (13) 

For completeness, the full expression for the energy per baryon is now displayed in terms of 
the individual nuclear, electronic, plus lattice contributions: 



e{A, Z] n) = J + xfVf [x^ + ylj - ln(xF + Vf) - Cej^Pp . (14) 

Note that given A, Z, and n = Atota\/V, the only unknown quantity in the above expression 
is the nuclear mass M{N, Z). While experimentally available for a large number of nuclei 
around the line of stability, nuclear masses near the drip line are unknown, thereby making 
the need for theoretical extrapolations unavoidable. As crustal properties become better 
determined, nuclear masses at the drip line will be strongly constrained. Alternatively, 
the advent of facilities capable of producing beams of rare isotopes to explore the limits of 
nuclear existence will place strong constraints on crustal properties. 

Having computed the energy per baryon of the system, we are now in a position to 
compute two additional thermodynamic properties that are essential for the understanding 
of both the structure and composition of the outer crust. These are the equation of state, 
namely, the relation between pressure and density, and the chemical potential. Recall that 
in modeling the outer crust the central assumption is that of thermal, hydrostatic, and 
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chemical equilibrium. Thus, complete equilibrium demands the equality of temperature, 
pressure, and chemical potential at each layer of the outer crust. 

At zero temperature and for a constant number of particles, the pressure of the system 
may be computed from the total energy of the system. As the individual nuclei do not 
contribute to the pressure, one must only compute the electronic and lattice contributions. 
In particular, the electronic contribution at zero temperature is given by 

fdEA xf fd^e\ ^ mf f ^ 3 r / n o\ , , j\ 

TZo-[-dv), = Y[drj~^' = ^^ V''' ' 8 r l^^ + y^) - ^"^"^ + ^^)J ) ■ 

(15) 

Similarly, the lattice contribution to the pressure is given by the following simple expression: 
In this manner the full {electronic plus lattice) contribution to the pressure may be written 



as 



(17) 



As alluded earlier, full equilibrium in the system is established by demanding that the 
temperature, pressure, and chemical potential — but not necessarily the baryon density — 
be continuous throughout the outer crust. As the temperature of the system is assumed to 
be equal to zero, the only remaining thermodynamic observable to calculate is the chemical 
potential. At zero temperature, the Gibbs free energy and the total energy of the system 
are related by a Legendre transform [G = E^PV\ That is, 

, G(A,Z-P) ,^ , P M(N,Z) Z 4:^ Z'^ 
f,{A, Z; P) = ' = e{A, Z;n) + -= ' + -/i^ - .Cej^PF , (18) 



where He = \/ v\e + "^e is the electronic chemical potential. Note that the chemical potential 
is a function of the pressure whereas the energy per baryon is a function of the baryon 
density. The transformation from one into the other is accomplished through Eq. (fT7|) . Also 
note that as hydrostatic and chemical equilibrium must be maintained throughout the star, 
the composition of the outer crust must be obtained by minimizing the Gibbs free energy 
per particle (i.e., /i) at constant pressure rather than by minimizing the energy per particle 
at constant baryon density. This procedure will be carried out in the next section. 



III. RESULTS 



In this section results will be presented for the structure and composition of the outer 
crust. The implementation of the ideas developed in the previous section will be carried 
out by using various models for nuclear masses. Two of these models are based on sophis- 
ticated microscopic/macroscopic models that yield root-mean-square (RMS) errors of only 
a fraction of an MeV when compared to large databases of available experimental nuclear 
masses ISj 



38| . These two models are the ones by Duflo and Zuker 39, 40, 4l[ and the finite 
range droplet model of MoUer, Nix, and collaborators 42|, |43|. The other two models are 
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based on accurately-calibrated microscopic approaches that employ a handful of empirical 
parameters to re pro duce the ground-state properties of finite nuclei and some collective ex- 
citations [H, While successful, the RMS errors of these two microscopic approaches 
are significantly larger than those obtained with the microscopic/macroscopic models. Yet 
one of the great advantages of the microscopic models is the ability to systematically study 
the impact of unknown physics on crustal properties. First and foremost, we are interested 
in understanding how models that are equally successful in describing available ground-state 
properties of finite nuclei differ in their predictions of exotic (neutron-rich) nuclei. 

A. Toy Model of the Outer Crust 

Although the structure and composition of the crust will be ultimately computed using 
sophisticated mass formulas, we start by introducing a "toy model" that while simple, cap- 
tures the essential physics of the outer crust, namely, a competition between an electronic 
density that drives the system towards more neutron-rich nuclei and a nuclear symmetry 
energy that opposes such a change. 

The toy model consists of the following two approximations. First, a simple liquid- 
drop model will be used to compute nuclear masses [see Eq. ([2])]. Second, the electronic 
contribution will be assumed to be that of an extremely relativistic [i.e., rrie/pFe^O) Fermi 
gas. While both of these approximations will be relaxed in the next section, we believe that 
the physical insights that one develops from this analytic treatment are valuable. 

The liquid-drop mass formula may be written in the absence of pairing correlations (and 
assuming Z{Z — 1)^Z'^) as follows: 

£n{x, y) = rupy + m„(l - y) - + — + a^x^y"^ + a.^{l - 2yf , (19) 

where x = A^/^ and y = Z/ A is the proton (or electron) fraction. The various empirical 
constants ( and Ca) represent volume, surface. Coulomb, and asymmetry contri- 

bution, respectively. Using a least-square fit to 2049 nuclei (available online at the UNEDF 
collaboration website http : // www . unedf . org/ ) one obtains the following values for the four 
empirical constants: 

= 15.71511 MeV, = 17.53638 MeV, = 0.71363 MeV, = 23.37837 MeV. (20) 

To understand the competition among the various terms — and to anticipate how this 
competition will be modified in the presence of a Fermi gas of electrons — we compute the 
optimal values of x and y using the simple liquid-drop formula by setting both derivatives 
equal to zero. That is, 

^ + 2a,xy'' = , (21a) 
Am + 2acX^y - Aa^{l - 2y) = , (21b) 
where we have defined Am = mn—mp. The above set of equations has the following simple 
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analytic solution: 



y = 7 — \ — ~ 7 — \ — • i^^u; 

\2aJ \AaJ 

The above solutions suggest the following physical interpretation. For a fixed proton fraction 
y = Z/A, the optimal value of x emerges from a competition between surface (which favors 
large x) and Coulomb contributions (which favors small x) . For the set of empirical constants 
given in Eq. fl2U]) . the relevant ratio is given by as/2ac~ 12.287. Conversely, if y4 = x^ is held 
fixed, then the optimal proton fraction y results from the competition between Coulomb 
and asymmetry contribution, with the former favoring y = and the latter y = l/2. If both 
equation are solved simultaneously, then one finds the most stable nucleus for this parameter 
set. One obtains xo = 3.906 and ?/o = 0.454, or equivalently: 

Ao = 59.598, Zq = 27.060, Nq = 32.538, (S/A)o = 8.784 MeV, mo = 930.195 MeV, (23) 



1/2 



(22a) 



with mo = (M/y4)o being the nuclear mass per nucleon. 

The second assumption defining the toy model is that of an extremely relativistic Fermi 
gas of electrons [i.e., p^e ^^g). In this limit one obtains simple expressions for the total 
energy per baryon, chemical potential, and pressure in terms of the adopted set of variables. 
That is, 

3 

£{x,y,p^) = en{x,y) + -y^^^Pp ~ Cex'^y^pp , (24a) 

4 

y, Pf) = £n{x, y) + y^'^PY - ^ Q x^v'^Vf > (24b) 
P(x, y, pp) = jy^^% - ^Ce x^y^pp . (24c) 



Before assessing the quantitative impact of the density dependent {i.e., electronic and 
lattice) contributions on the semi-empirical mass formula, a few comments are in order. 
First, at the predicted neutron-drip density of about 4xl0^^g/cm^, the Fermi momentum is 
approximately equal to pp'^^^40 MeV. This suggests a large electronic contribution at those 
densities of about e™^"" ~ 30 y'^^^ MeV. As the nuclear contribution is independent of density, 
the electrons will drive the system to small values of y. Second, the lattice contribution 
(perhaps not surprisingly) has the same dependence on x and y as the Coulomb contribution 
to the semi-empirical mass formula. Indeed, the full impact of the lattice contribution can be 
included through a redefinition, albeit a density dependent one, of the Coulomb coefficient. 
That is, ttc ^cicipp) = (ac — C^pp). As the optimal value of the proton fraction y (for fixed 
x) emerges from a competition between Coulomb and asymmetry terms [see Eq. (]22b|) ]. the 
lattice contribution drives the system towards the symmetric y = l/2 limit. Yet the lattice 
contribution is marginal. Indeed, even at neutron-drip densities its contribution provides a 
meager, although by no means negligible, correction to the dominant Coulomb term. The 
above two facts summarize the main structure of the outer crust, namely, a nuclear lattice 
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embedded in an electron gas that is responsible for driving the system towards progressively 
more neutron-rich nuclei. Thus the outer crust represents a unique laboratory for the study 
of neutron-rich nuclei in the Z^20—50 region. As such, it nicely complements Rare-Isotope 
Facilities worldwide that aim to provide a detailed map of the nuclear landscape. 

Incorporating electronic and lattice contributions to the semi-empirical mass formula 
yields the following expression for the total energy per nucleon of the system: 

e{x, y, pp) = rripy + m„(l - y) - + — + a^x^y^ + a^{l- 2yf + -^y^'^Pp - Ci x^y^pp . (25) 

X 4 

As done before for the pure nuclear contribution, the optimal values of x and y — at 
fixed density — may be obtained by setting both derivatives equal to zero. That is. 



de\ a. 
dx 

de 



dx j X 



2 + 1a,xy' = , (26a) 



= -Am + 2a,x'y - Aa^{l - 2y) + = , (26b) 



dy 

where a "renormalized" Coulomb coefficient has been defined as 

«c(Pf) = (etc - Cepp) . (27) 



1. First- order Solution 



Before providing exact solutions to Eqs. we compute approximate solutions that 

are accurate to first order in pp. In addition of being analytic, these closed-form expressions 
provide valuable insights into the composition of the outer crust. The first-order solutions 
are obtained by incorporating the density dependence in the following form: 



x{Py) = a;o(l + and y{pp) = yo{l + v) 



(28) 



where both ^ and rj represent small {i.e., first-order in pp) deviations from the zero-density 
results. Substituting the above equations into Eqs. fl2Bl) yields the first-order solutions. One 
obtains, 



x{pp) 



Xq 



yo 



1 + 



(Ci - 1)C, + 2C2 

3Ci - 1 
3C2 -CA pf 
3Ci — 1 J 



Pf 



(3.90610 + 0.03023pf) 
(0.45405 - 0.00419pf) . 



(29a) 
(29b) 



Note that in the above expressions the Fermi momentum should be given in MeV. Moreover, 
for simplicity the following two dimensionless quantities were introduced: 



2 

XqOiq 



8.58843 and 



1 



0.05547 . 



(30) 



The first-order equations [Eqs. fl29l) ] — while not necessarily quantitatively accurate — 
provide useful insights into how the composition of the outer crust evolves with density. As 
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previously suggested, the proton fraction y decreases with density in an effort to minimize 
the "repulsive" electronic contribution. Indeed, to an excellent approximation Eq. (129bp 
may be written in the following simple form: 

1/(Pf) =yo-^ = (0.45405 - 0.00411pp) . (31) 

As indicated in Eq. (122b|) . the optimal value of emerges from a competition between 
Coulomb and asymmetry terms, with the former driving y^ towards zero and the latter 
towards one half. The above equation indicates that the evolution of y with density is 
controlled by the dimensionless ratio of p-p^/'^a, suggesting that the larger the value of the 
asymmetry energy, the slower the evolution away from ?/o; that is, the more symmetric 
the nucleus will remain. Moreover, as the denominator in Eq. (13 ip [Soa ~ 100 MeV] is 
significantly larger than the electronic Fermi momentum over the entire region of interest, 
the first-order approximation is expected to be fairly accurate over the entire outer crust. 
Indeed, assuming a realistic value for the drip density of pdrip = 4xl0^^ g/cm yields a proton 
fraction of |/jjj.ip = 0.298. This represents a 2% deviation from the value of y{}^^¥^i) = 0.305 
for the conventionally accepted drip nucleus ^^^Kr. 




10 20 30 40 



Pp=(37t^n)^^^ (MeV) 

FIG. 2: (color online) Baryon number x = A^^^ (lower panel) and proton fraction y = Z/A (upper 
panel) are displayed as a function of the Fermi momentum = (37r^n)^/^. The black solid lines 
represent the exact solution to the toy-model problem given in Eqs. (I26p . while the red dashed 
lines display the corresponding solution in the C^ = (no lattice) limit [see Eqs. (j32p ]. Finally, the 
low-density solution [Eqs. ([29P ] is displayed by the blue dot-dashed lines. 
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2. Exact Solution 



The exact solution to the toy-model problem requires (for a fixed value of p-p) to find the 
simultaneous roots of Eqs. fl26ap and fl26bp . While numerically simple, the exact solution 
can not be displayed in closed form. Yet the exact solution differs only slightly from the 
Ci = solution — which has an analytic, albeit a bit unorthodox, solution. The closed-form 
solution for the Ci = case may be obtained by simply re- writing Eqs. fl2Bl) . That is, 

pM = + . (32b) 

The above set of equations suggest that rather than looking for a solution of x and y as 
a function of pp, one should "solve" for x and function of y, with the maximum 

value of y given by ?/max = Z/o = 0.45405 and the minimum value of y given by the condition 

/^(l/min) 

In Fig. [2] the baryon number x = A^^^ and proton fraction y = Z/A are displayed as a 
function of the Fermi momentum pp = (Svr^n)^/'^. The black solid line display the exact 
numerical solution to the toy-model problem [see Eqs. (126!) ]. In this simple model, the drip 
line density is predicted to be at pdrip = 4xl0^^ g/cm^ with the drip-line nucleus being ^^'^Cd 
{i.e., Z = 48 and N = 106). The solution obtained by ignoring the lattice contribution is 
displayed by the red dashed line. Because the lattice contribution to the chemical potential is 
negative, the Ce = solution reaches the drip line faster, i.e., at a lower density. Moreover, 
as the lattice contribution "renormalizes" the Coulomb term in the semi-empirical mass 
formula (or equivalently, enhances the role of the symmetry energy) the Ce = solution 
predicts a lower proton fraction than the exact solution. Finally, the dot-dashed blue line 
displays the solution correct to first-order in p^. In the particular case of the proton fraction 
y, the approximate linear solution y = yo —ppe/^^a [Eq- fl3T]) ] reproduces fairly accurately 
the behavior of the exact solution. 

The equation of state {i.e., pressure vs density) predicted by the toy model is displayed 
in Fig. [31 As the lattice provides a negative contribution to the pressure [Eq. fl24cl) ]. the 
equation of state for the Ce = case is slightly stiffer than the exact one. The first-order 
solution in pp provides a quantitatively accurate description of the equation of state up to 
fairly large values of the density. Note that the first-order approximation to the pressure is 
defined as follows: 

-jy^^^ - x^y'^ ^ (0.08367 - 0.00106^^) . (33) 



npp 4 3 



B. Realistic Models of the Outer Crust 



In this section we employ realistic nuclear mass models to compute the structure and 
composition of the outer crust. Two of the models 



40, 41, 43 are based on so- 



phisticated mass formulas that have been calibrated to thousands of available experimen- 
tal masses throughout the periodic table [13, 3^- The other two models are based on 
accurately-calibrated microscopic approaches that employ a handful of empirical parame- 
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p^=(3Knf^ (MeV) 

FIG. 3: (color online) Pressure as a function of the Fermi momentum pp = (Svr^n)^/'^. The black 
solid line represent the exact solution to the toy-model problem, the red dashed line displays the 
corresponding = (no lattice) solution, and the low-density solution is displayed by the blue 
dot-dashed line. 



ters to reproduce the ground-state properties of finite nuclei and some nuclear collective 



excitations ^2|, |13|, [14 



Whereas the microscopic models have yet to attained the level of precision displayed by 
the microscopic/macroscopic ones, they are valuable in elucidating various details of the 
underlying physics. For example, in the previous section we established the critical role 
played by the symmetry energy in the evolution of the proton fraction with density [see 
Eq. (l3T!l ]. However, it is unknown how the symmetry energy coefficient changes as nuclei 
move towards the drip line. Presumably, the development of a significant neutron skin 
makes these nuclei (on average) more dilute than their stable counterparts. If so, one needs 
to extrapolate the symmetry energy to lower densities, a procedure that is highly uncertain 
because of our poor knowledge of the slope of the symmetry energy. To illustrate this 
uncertainty, the equation of state of pure neutron matter predicted by NL3 (green dashed 
line) and FSUGold (blue solid line) is displayed in Fig. |H For comparison, we also show the 
predictions from the microscopic model of Friedman and Pandharipande based on realistic 
two-body interactions ^] (purple upside-down triangles) and the model-independent result 
based on the physics of resonant Fermi gases by Schwenk and Pethick (red hatched 
region). Note that to a very good approximation, the equation of state of pure neutron 
matter equals that of symmetric nuclear matter plus the symmetry energy. The differences 
between NL3 and FSUGold displayed in Fig. H] are all due to the large uncertainties in the 
symmetry energy. In particular, as NL3 predicts a stiffer equation of state than FSUGold, 
namely, one whose energy increases faster with density at high densities, the symmetry 
energy of NL3 is lower than that of FSUGold at sub-saturation densities. Thus, FSUGold 
has been shown to reach the neutron-drip lines earlier than NL3 [i^. By the same token, 
NL3 should predict a sequence of more neutron-rich nuclei (lower y) in the outer crust than 
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FIG. 4: (color online) Energy per particle for pure neutron as a function of the neutron Fermi mo- 
mentum. Shown are the microscopic model of Friedman and Pandharipande [4j] (purple triangles) 
and the model-independent result based on the physics of resonant Fermi gases by Schwenk and 
Pethick 0] (red region) . Also shown are the predictions from the accurately calibrated NL3 
(green dashed line) and FSUGold [ij] (blue line) parameter sets. 
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FSUGold. 

Shown in the left-hand panel of Fig. [5] is the proton fraction predicted by the two mi- 
croscopic models; FSUGold (blue solid line) and NL3 (green dashed line). Also shown is 
the simple prediction obtained from the liquid-drop formula [Eq. (!3T!) ]. The proton fraction 
predicted with the FSUGold parameter set is consistently higher than for the NL3 set. This 
is a reflection of the stiffer penalty imposed on the FSUGold set for departing from the 
symmetric {N = Z) limit. The right-hand panel shows the corresponding behavior for the 
case of the microscopic/macroscopic models of MoUer-Nix (red solid line) and Duflo-Zuker 
(purple dashed line). Differences among these two models are small. 

Similar trends may be observed in Figs. Eland [7] where the composition of the outer crust 
is displayed as a function of density. As the system makes a rapid jump in neutron number 
(say to magic number N = 50) the proton number jumps with it. Along this neutron plateau, 
the proton fraction decreases systematically with increasing density in an effort to reduce the 
electronic contribution to the chemical potential. Eventually, the neutron-proton mismatch 
is too large and the jump to the next neutron plateau ensues; a jump that is driven by 
the symmetry energy. Clearly, the larger the symmetry energy at low densities, the smaller 
the neutron-proton mismatch and the early the jump to the next neutron plateau. These 
features are clearly displayed in Fig. [H] as one contrasts the behavior of FSUGold to that 
of NL3. In contrast, few differences are noticeable in Fig. [7] when comparing the model of 
Moller-Nix to that of Dufio-Zuker. 

We conclude this section by displaying in Fig. [8] equation-of-state (pressure vs density) 
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MoUer-Nix 
Duflo-Zuker 

y=yo-PFe/8a, 
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Pp=(37l^n)''^ (MeV) Pp=(37i^n)"^ (MeV) 

FIG. 5: (color online) Shown in the left-hand panel is the proton fraction predicted by the accurately 
calibrated FSUGold (blue solid line) and NL3 (green dashed line) parameter sets. Also shown is 
the simple liquid-drop formula given in Eq. (|3ip . The right-hand panel displays the corresponding 
proton fraction as predicted by the mass formulas from Moller-Nix (red solid line) and Duflo-Zuker 
(purple dashed line). 



predictions for the outer crust of a neutron star. The left-hand panel shows results from 
calculations using the FSUGold (blue solid line) and NL3 (green dashed line) parameter 
sets. Although barely visible, the density shows discontinuities at those places where the 
composition changes abruptly. It is also noted that the FSUGold parametrization predicts 
a pressure that rises slightly faster with density than NL3. For the NL3 set, the symmetry 
energy admits lower values of the proton/electron fraction y which, in turn, lowers the 
pressure of the system. Lower values of y also yield lower values of the chemical potential, 
thereby delaying the arrival to the neutron-drip line. Indeed, whereas FSUGold predicts a 
drip-line density of pdrip = 4.17xl0^^g/cm'^, with NL3 the transition is delayed by about 8%, 
or until pdrip = 4.49 x lO^^g/cm'^. A similar plot is shown for the microscopic/macroscopic 
models of Moller-Nix (red solid line) and Duflo-Zuker (purple dashed line). Differences 
among these two models are barely noticeable. Indeed, drip-line densities in both models 
are predicted at about pdrip = 4.3 x lO^^g/cm^. Model predictions for various observables at 
the base of the outer crust {i.e., in the drip-line region) are listed in Table [B 



IV. CONCLUSIONS 

Following the seminal work by Baym, Pethick, and Sutherland, as well as the more recent 
comprehensive work by Ruester, Hempel, and Schaffner-Bielich, we studied the composition 
and equation of state of the outer crust of non-accreting neutron stars. The central focus of 
our study was the sensitivity of crustal properties to the density dependence of the symmetry 
energy. To do so, four different models were adopted. Two of these models, Moller-Nix 
and Duflo-Zuker, are based on a combined microscopic/macroscopic approach and yield 
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FIG. 6: (color online) Composition of the outer crust of a neutron star as predicted by the accurately 
calibrated FSUGold (upper panel) and NL3 (lower panel) parameter sets. Protons are displayed 
with the (lower) blue line while neutrons with the (upper) green line. 



Model 


p(10"g/cm^) 


n(10-^fm-^) 


P(10-'^MeV/fm^) 


/ie(MeV) 


Element 


B/A{Me\) 


Moller-Nix 


4.34 


2.60 


4.93 


26.22 


iiSKr 


7.21 


Duflo-Zuker 


4.32 


2.58 


4.89 


26.17 


iiSKr 


7.19 


FSUGold 


4.17 


2.50 


4.68 


25.88 


iiSKr 


7.11 


NL3 


4.49 


2.69 


5.06 


26.39 


i20Kr 


7.13 



TABLE I: Equation-of-state observables (mass density, baryon density, pressure, and electronic 
chemical potential) and composition (nucleus and binding-energy per nucleon) at the base of the 
outer crust. 



the most accurate nuclear masses available in the literature. The other two models, NL3 
and FSUGold, are of a purely microscopic nature and based on a relativistic mean-field 
approach. Although the former are significantly more accurate than the latter, microscopic 
models have the advantage of making definite predictions on how the symmetry energy 
changes with density (see Fig. H]). One can then study the impact of various features of the 
symmetry energy — such as its slope — on crust al properties. 

The composition and equation of state of the outer crust emerge from a competition 
among three relatively simple contributions to the total energy (or chemical potential) of 
the system: nuclear, electronic, and lattice. The nuclear contribution appears exclusively 
in the form of nuclear masses and is independent of the baryon density. The electronic 
contribution is modeled after a zero-temperature free Fermi gas and dominates the behavior 
of the system with baryon density. Finally, the (body-centered-cubic) lattice contribution is 
also density dependent and provides a relatively modest correction (of no more than 10%) to 
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FIG. 7: (color online) Composition of the outer crust of a neutron star as predicted using the the 
mass formulas of Moller-Nix (upper panel) and Duflo-Zuker (lower panel). Protons are displayed 
with the blue (lower) line while neutrons with the green (upper) line. 



the energy of the system. The alluded competition is then primarily driven by an electronic 
term that favors a small electron fraction (to reduce the electronic Fermi energy) and a 
nuclear symmetry energy that opposes such a shift towards progressively more neutron- 
rich nuclei. To motivate the discussion and to highlight this competition, we implemented 
a "toy model" of the outer crust by using a simple semi-empirical ( ^^Bethe-Weizsackei") 
nuclear mass formula. Volume, surface. Coulomb, and asymmetry terms were extracted 



from a least-square fit to 2049 nuclei (see http : / /www . unedf . org/ ) . The advantage of such 



a simple model is that useful insights emerge from the analytic structure of our results. 
Indeed, a particularly transparent result that illustrates nicely the competition between the 
electronic contribution and the nuclear symmetry energy was obtained, namely, 

= 2/0 - + 0(PFe) , (34) 

where yo is the zero-density proton fraction, p^e is the electronic Fermi momentum and 
Oa is the symmetry energy coefficient. While illuminating, this (first-order) result is also 
surprisingly accurate, as the electronic Fermi momentum at the base of the outer crust is 
very close in value to the symmetry energy coefficient (ppg ~ 26 MeV vs aa ~ 23 MeV). In 
particular, the toy model predicts a value for the electron fraction at the base of the crust 
that differs by only a few percent from that of the drip-line nucleus ^^^Kr. 

What is unknown, however, is how the symmetry energy coefficient is modified as 
nuclei move away from the line of stability. Presumably, the symmetry energy is reduced 
in neutron-drip nuclei due to the development of a dilute neutron skin. To investigate the 
sensitivity of the outer crust to the density dependence of the symmetry energy we employed 
two relativistic mean-field models (NL3 and FSUGold) that while accurately calibrated, 
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FIG. 8: (color online) Shown in the left-hand panel is the zero-temperature equation of state 
(pressure vs density) predicted by the accurately calibrated FSUGold (blue solid line) and NL3 
(green dashed line) parameter sets. Also shown is the prediction from the simple liquid-drop 
formula. The right-hand panel displays the corresponding expression as predicted by the mass 
formulas from Moller-Nix (red solid line) and Duflo-Zuker (purple dashed line). 



predict a significantly different density dependence for the symmetry energy. In particular, 
NL3 predicts a smaller symmetry energy than FSUGold at the (small) densities of relevance 
to the outer crust (see Fig. H]). One of the main goals of the present manuscript was to 
document how such differences impact the composition of the outer crust. One noticed, 
quite generally, that as the density increases along a fixed neutron-number plateau (say at 
magic number N = 50) the proton fraction decreased systematically in an effort to reduce the 
electronic contribution to the chemical potential. Eventually, however, the proton fraction 
becomes too low and the symmetry energy drives the system into the next plateau (say at 
magic number A^ = 82). How low can the proton fraction get is then a question that must be 
answered by the symmetry energy. Indeed, whereas NL3 predicts the formation of 28Ni5o, 
FSUGold (having a larger symmetry energy) leaves the = 50 plateau with the formation 
of 32Ge5Q (or four protons earlier). This result may be stated in the form of a correlation 
between the neutron radius of ^°^Pb and the composition of the outer crust: the larger the 
neutron skin o/^°^Pb, the more exotic the composition of the outer crust. Finally, and as it 
was done in Ref. [29|, we have computed crustal properties using two of the most accurate 
tables of nuclear masses available today, namely, those of Moller-Nix and Duflo-Zuker. Our 



results using the model of Moller and Nix agree well with those published in Ref. [29 
These results are practically indistinguishable from the ones obtained using the Dufio-Zuker 
nuclear mass table; a table that includes 9210 nuclei! 

In summary, we have used realistic nuclear mass tables to elucidate the role of the sym- 
metry energy on the structure and composition of the outer crust of neutron stars. Recent 
observations of crustal modes in magnetars are likely to provide stringent limits on the 
equation of state of neutron-rich matter. As the field of nuclear astrophysics continues to 
advance — with the commission of both radioactive beam facilities as well as ground- and 
space-based telescopes — we enter a new era that promises great hope in the determination 
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of the nuclear matter equation of state. 
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